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Abstract 

We develop the quantum theory of an open-cavity laser assuming that only two modes compete 
for gain. We show that the modes interact to build up a composite mode that becomes the lasing 
mode when pumping exceeds a threshold. This composite mode exhibits all the features of a 
typical laser mode, whereas its precise behavior depends explicitly on the openness of the cavity. 
We approach the problem by using the density-matrix formalism and derive the master equation 
for the light field. Our results are of particular interest in the context of random laser systems. 

PACS numbers: 42.55.Ah, 42.55. Sa, 42.55.Zz 
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I. INTRODUCTION 



Small size, complex structure and extreme openness or complete absence of the cavity 
are characteristic features of a number^f "exotic" laser systems that have attracted the 

Examples of lasers that fall into this category 



attention of physicists in recent years 111. 

n n 

are chaotic microcavity [2|, |3| and random 



J-[6j lasers. These systems are quite different 



from "traditional" cavity lasers composed of an amplifying medium in a high quality-factor 
cavity 0, [sl . From the theoretical point of view, the very strong coupling to the external 
world requires a special treatment that is different from what can be found in standard laser 
textbooks [3, Isl. A semiclassical model of lasing in open complex or random media was 
developed by Tiireci et al. This theory was successfully applied to understand lasing 



in random media 



10| . To develop the quantum theory of an open-cavity laser, one starts 



by facing the problem of quantization of the electromagnetic field in a space that cannot 
be separated into "system" and "bath" parts unambiguously. This problem was solved by 



Hackenbroich et al. 



11| who also put forward Langevin and master equations to describe the 



dynamics of modes in open resonators 



121 ]. Hackenbroich also derived Heisenberg-Langevin 



equations for an open-cavity laser, which, however, he analyzed only in the semiclassical 
approximation 



13| . A related problem of light emission by an atom in a lossy cavity was 
also considered by Di Fidio et al. [l^ who, however, considered a single-mode field and 
hence didn't discuss additional features that arise from the coupling between modes due to 
the openness of the system. 

In the present paper we use a combination of the quantization procedure of Refs. 11 1 
with the standard density operator approach to develop the full quantum theory of an 
open-cavity laser. We compute and analyze the lasing threshold, the photon statistics (both 
below and above the threshold), as well as the emission linewidth of a laser that has no 
well-defined cavity, assuming that only two modes compete for gain. This simple two-mode 
model allows us to capture some of the essential features of cooperative mode dynamics that 
seems to determine the behavior of the system. We compare our results with those known 
from the standard laser theory js] and highlight common features as well as important 
differences. It is worthwhile to note that a different master equation for a random laser 
was previously proposed by Florescu and John 15|. These authors considered the random 



laser as a collection of low quality-factor cavities, coupled by random photon diffusion. In 
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contrast to this work, our approach has the advantage of not relying on any particular model 
of wave transport, as well as being based on a well-defined quantization procedure for the 
electromagnetic field and a fully quantum model for the atoms providing amplification. 



II. MASTER EQUATION FOR THE REDUCED DENSITY OPERATOR OF THE 
ELECTROMAGNETIC FIELD 

Let us start by considering an ensemble of two-level atoms interacting with the electro- 
magnetic field. We divide the modes of the electromagnetic field into those that belong to 
the system "atoms + field" (A + F) and those that constitute the "bath". In the density 
operator approach, the system A + F is described by the density operator p{t). The reduced 
density operator prit) describing the electromagnetic field (F) is obtained by tracing over 
the atomic (A) degrees of freedom: ppit) = Tr^/3(t). The dynamics of the laser is due to 
the competition between gain (due to the interaction of the field with atoms) and loss (due 
to the coupling of the system A -|- F to the bath), which can be expressed in the form of 
the following master equation: 

= n^^'-^'^pF + L(^°^^)p^. (1) 

Here the super-operators 

^(gain) ^(loss) describe the gain and the loss, respectively. 
Equation ([1]) is quite formal and can be written for any quantum system interacting with 
environment. Let us now give expressions for the two terms on the r.h.s. of Eq. ([1]) in the 
case of an open-cavity laser. A way to deal with the second term L^^^^^'^pp was proposed 



in a series of papers by Hackenbroich et al. 11 



12| . The idea is to separate the physical 



space into two subspaces and to quantize the field in terms of the modes a and h of these 
subspaces. In the context of the laser system considered here, it is natural to choose the first 
subspace such that it contains all the atoms and has a finite volume. The discrete modes 
of the first subspace will constitute our sub-system F, whereas the modes of the second 
subspace will make up the bath. An equation for the density matrix of the sub-system F is 
derived by tracing over the degrees of freedom corresponding to the modes that belong to 
the bath. This yields Q 

^(loss)^^ = ^7aa' {2ax'pF(^\ - pFa{ax' - a\ax'pF^ ■ (2) 

A,A' 



Here d\ and a\ are annihilation and creation operators corresponding to the modes of the 
sub-system F, and the coefficients 7aa' depend on the precise geometry of th e sy stem. These 



coefficients were calculated for a number of particular open cavities in Ref. [16|] but may be 
difficult to obtain in the general case. In a random laser system, they may be treated as 
random variables [l^. 

The essential difference between Eq. ([2]) and the analogous equation of the standard laser 
theory is that the damping matrix 7 is not diagonal. This shows that the openness of 
the system not only leads to losses described by the diagonal elements of 7, but also induces 
coupling between different modes. The strength of the coupling is given by the off-diagonal 
elements of the damping matrix 7. 

Let us now turn to the ffist term in Eq. ([1]). It is not specific for the op en- cavity laser, so 
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. As 



that we will follow standard approaches to derive an explicit expression for it 
the ffist step, we consider the Jaynes-Cummings Hamiltonian for an atom interacting with 
the electromagnetic field (we set h = 1 in the following) [l^ : 

H = + ^ uJxa{ax + ^{gxcr^ax + h.c). (3) 

A A 

Here Ua is the frequency of the atomic transition, = \e){g\ and az = |e)(e| — \g){g\ are the 
atomic raising and inversion operators, respectively, the states \g) and |e) are the ground 
and excited states of the two-level atom, ux are the frequencies of the modes of the field, and 
the coefficients gx describe the coupling between the atom and the mode A of the field. It 
is convenient to introduce a reference frequency u and the detuning parameters 6 = Ua — 
and Ax = ux — 00 to write 

Co S 
H = —Gz + ^ ^ a\ax + -0-2 + ^ Axa\ax + ^(^'AO-'^flA + h.c.) 

A A A 

= ^0 + V^, (4) 

where -ffo = ^o'z/'^ + ^ XIa ^a^a and V = 6az/2 + J2x ^A^j^aA + J^xid^^'^^^^ + h.c). Because 
Ho and V commute, [Hq, V] = 0, we will work in the interaction picture where the dynamics 
of the system is governed by V. In this picture, the time evolution of the density operator 
is given by the evolution operator U{t) = exp[—iVt]: 

pit) = u{t)mu\t), (5) 
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where 

We will now restrict our consideration to the situations when the frequencies ux of the modes 
are close to oj, such that ^ gx, S, and will proceed by setting Ax = 0: 

V - ] . (T) 

Here gA = Ea^'a^a and g = (EaS'a)^^^- "^^^ newly defined operator A obeys the standard 
bosonic commutation relation [A, A^] = 1. After some algebra, with the help of the operators 
= giAA'^ + {6/2gy]^/'^ and = g[A^A + {6/2gYY^'^, the evolution operator reads 

^^^^ ^ j cos[^t] -i(5/2)sin[(^t]/(^ -ism[ipt]/ipA \ 

sin[(pt]/(p cos[</)t] + i (5/2) sin[0t]/0 / 

We now assume that at the initial time t = 0, /)(0) = Pi?(0) ® Pa(0), with being the 
density operator of the atom, and that the atom is in its upper state: Pa(0) = |2)(2| = 
{(7z + l)/2. The density matrix of the full system "atom + field" at t = is then 



p(o)=r'-'"K (9) 




Considering ^± = cos[0t] ± i (5/2) sm[0t]/(p, equations (jS]) and ([HD allow us to compute the 
density operator at arbitrary time as 



^ <I_Pf(0)<I+ ig^_pF{O)sm[0t] /0A ^ 



Pit) 



(10) 



^— i^'A"!' sin [0t] /0Pp{O)^+ g'^A^ sin [0t] /0pp{O) sin [(^t] /(^A^ 
Finally, the reduced density operator ppit) = Ti Apit) is 

prit) = ^^pf{0)^+ + g^A^ sm[^t] /^pF{0)sm[(pt] /(pA 

= A(t)PF(0), (11) 

where we defined a super-operator A(t). 

Equation (fTTj) yields the evolution of the reduced density operator of the electromagnetic 
field interacting with a two-level atom which is initially in the excited state. This is clearly 



insufficient to describe laser emission. The first lacking ingredient stems from the fact that 
we want to describe an ensemble of many atoms, not just a single atom. To make use of 
Eq. f llip in the case of many atoms, we assume that (i) the atoms interact with the field 
one after another, in sequence, and not all at a time, and (ii) the time of interaction of a 
given atom with the field r is much shorter than the typical time t at which the evolution 
of the field is calculated. The density matrix of the field after a time t ^ t during which 
the field interacted with k atoms will be then equal to pp\t) = A'^{t)pp{0) js], [isl. We 
will now introduce the second important ingredient of the laser system — the pump. To 
model the pump in the framework of our two-level atom model, we assume that atoms 
are transferred to the excited state at a rate r by some external mechanism (for example, 
via additional atomic levels that are not included in our model explicitly) and that the 
probability for k atoms to get excited during a time At is P{k) = CKkp'^i^ — p)^~^i where 
Cxk = K\/k\{K — k)\, p is the probability for a given atom to be in the excited state, and 
K is the total number of atoms that can potentially participate in the lasing process (i.e., 
< A; < K). The average number of excited atoms is {k) = pK = rAt. The parameter p 
describes statistics of pumping, with the limit p — > (that we will consider from here on) 
correspondin g to random pumping and the limit p — ?■ 1 corresponding to a uniform (regular) 



pumpmg 
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20|. 



The density operator averaged over k is 
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1^ 



^ K 



p^(t) = J]P(A:)p?)(t) = {l+p[A(r)-l]} p^(0). (12) 

fc=0 

To obtain a dynamic equation for pF{t), we take the derivative of Eq. (fT2l) with respect 
to time and expand the result in series in p{A — 1): 



Mt) = -ln|l+p[A(r)-l]}p^(t) 



p 

^ r[A(r) - l]ppit) - |[A(r) - l^Mt). (13) 

Finally, we now take into account the fact that the time of interaction of a given atom 
with the field r is, in fact, a random variable, r is finite due to the possible decay of the 
excited state without coupling to the modes of the electromagnetic field that make part of 
our sub-system F. This decay may be due, for example, to transitions involving additional 
atomic levels (with or without emission of a photon), not included in our model. With F 
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being the rate of such transitions, the statistical distribution of r is P{t) = T exp(— Fr). By 
averaging Eq. f fT3|) over this distribution, we obtain 

L^'^'^'^PF = r drV eM-^r) {[A(r) - 1] - |[A(r) - l]^}p^(t). (14) 

Equations ([2]) and f|T4|) provide exphcit expressions for the two terms on the r.h.s. of Eq. 
(II]). It is worthwhile to note that Eq. ([2]) was derived in the Schrodinger picture, whereas 
Eq. ( 1141) — in the interaction picture. In the interaction picture, the general form of Eq. ([2]) 
remains unchanged, except for the damping matrix that has to be transformed accordingly. 
In the present paper we will not use any particular model for this matrix but will rather 
treat it as a free parameter, having in mind that in a random laser, for example, it is a 
random matrix. 



III. TWO-MODE MODEL 

Under general conditions, many modes may coexist and compete for gain in an open- 
cavity laser. The off-diagonal elements of the damping matrix 7 couple the modes and 
make the single-mode regime hardly realizable. It is worthwhile to note that this coupling 
is different from the coupling via interaction with the atomic subsystem which is always 
present and hardly depends on the type of the cavity under consideration. To analyze the 
interaction between modes, we consider a model in which only two modes are taken into 
account. This simple situation often allows for important insights into the dynamics of laser 
systems, as, for example, it was the case for the quantum-beat or correlated-emission lasers 



2l|. Our model differs from the previously considered two- mode models (see, e.g., Ref. 21 1 



but also Refs. 



22| and |23|) by the mode coupling through both the atomic subsystem (i.e., 



the coupling due to the fact that the modes interact with the same atomic transition' 
the common bath, whereas only the first type of coupling was considered in Refs. 21 



and 



23|. 



Thus, having a common bath for the two modes is essential in our model. The strength of 
the additional coupling between the modes is given by the off-diagonal elements 712 and 721 
of the damping matrix 7. To put accent on this new element of the model, we will focus on 
the dependence of our results on these off-diagonal elements in what follows. To proceed. 
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we rewrite the master equation ([T]) assuming the hmit of p — > 0: 

POO |- -, 

p = r drV exp{—TT) ^^p{t)^+ + g'^ce'^ sin [(pr] /if p{t) sin [(pr] /ipa — p{t) 
Jo 

+ (15) 

where a = {giai + g2a'2)/g and cp = g[aa^ + {5/2g)'^Y^'^. To hghten the notation, we drop 
the subscript "F" of the reduced density operator and write pp = p. The hmit p 
corresponds to the Poissonian distribution of the number of excited atoms, P{k), and hence 
to the reahstic case of random pumping by an external source. The loss term in Eq. ( IT5|) 
follows from Eq. 

l^(ioss)^ = 7ii(2aipa| — pd{ai — d\dip) + 7i2(2a2paJ — pd\a2 — a[d2p) 

+ 72i(2aipd2 - pdlai - a\aip) + 722(202/002 - pa2a2 - d\a2p). (16) 



To solve the master equation f[T5|) . we will work with composite modes a defined above and 
13 = {g2ai — gia2)/g, with the commutation relations: = = 1 and [a, = 

[a, (3] = (here we assume that gi^2 are real numbers). The composite modes a and /3 are 
not due to the phase locking phenomenon but are simply linear superpositions of the bare 
modes 1 and 2. Using the modes a and /3 instead of 1 and 2 simplifies further analysis 
because only the mode a will be actually excited in the lasing process, as we will see from 
the following. 

To simplify analytical calculations, we will assume 721= 712 from here on. This corre- 
sponds, for example, to open cavities considered in Ref. 16|]. The density operator can be 
represented in the basis of Fock states |no,,n^) as 

P= ^ Pn^,np;m^,mp\na,np){mo,,mp\. (17) 
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Equations ( !T5|) and (fT6|) yield an equation for the density matrix Pna,np;mo,,mf,'- 

' A{na + + 2)/2 + Zy45(ra„ - mQ)/2 + B{na - maY/8 
_ 1 + P + {B/2A){n^ + m„ + 2) + (5/4A)2(n, - m,)^ 

+ Ci(?7,a + mQ,)/2 + C2(r;,/3 + m/3)/2 

+ C'lVl'^a + + l)Pn<,+l,n^;m<,+l,m^ + C2\J ijlp + l)(m/3 + l)pn<,,n^+l; mc,m^+l 

Wq — l,n^+l; ma,mp ) 

(18) 

where we defined A = 2r(^7/r)2, 5 = 4((7/r)2A, 5 = 5/r, Ci = 2g-\^^^gl+2^^29ig2+l22gl), 
C2 = 2g-^{-fng2- '^1129192 + 12291) and C3 = fi'"^[(7ii -722)fi'ifi'2 + 27i2(fi'i-fi'?)]. In order to 
facihtate the comparison with the standard laser theory, we defined the coefficients A and B 
in the same way as in the book js] (p. 333). When 712 = and 711 = 722, the coefficients Ci 
and C2 reduce to the coefficient C of this book (p. 255) and C3 vanishes. The coefficient C3 
vanishes also for 711 = 722 and gi = 92- This case is somewhat special because it corresponds 
to two modes with the same losses and the same coupling with atoms. This leads to C3 = 0, 
but still Ci 7^ C2 for 712 7^ and hence the problem does not reduce to the case of a high 
quality-factor cavity. High quality-factor cavities are described by Eq. ( fT8l) with 712 — )■ 0. 
The coupling between the bare modes 1 and 2 then arises uniquely from their interaction 
with the same atomic transition (and not with the common bath) and the composite modes 
a and /3 still provide a useful basis for the description of lasing (see, e.g., Ref. |2l[|). 



For the diagonal elements of the density matrix Pna,ng;na,ng = Priori p and up to the second 
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order in B /A = A^g/T)"^ ^ 1 we obtain 



Pna,np 



AUn 



2C|(n„ + l)np 2Cln^{np + 1) 



A{n^ + 1) 



+ CiUa + C2ni3 



.C|K + l)(r2^ + l) 

Pna,np + Pna+l,nj3+l 



na + l,na + l 



+ Pna+l,n0~l + PnQ"-l,n^+l 



+ 
+ 



4C|(n„ + l)(n;3 + l) 4C|(n„ + l)n^ 



CiK + 1) 



C2{np + 1) 



4C|(n„ + l)(r2^ + l) 4C3V(r^/3 + l) 



na,np+l 



Pna+l,nis 
Pna,np+l-i 



(19) 



where = M„,,,„^ + (M/2)2[l + 5^ + + 1/2) + (5/4A)2]-2m-J^„^ and 

M^^^n, = [A{n^ + 1/2) + fi/4][l + 52 + (5/A)K + 1/2) + (5/4A)2]"i + CiK " 1/2) + 
C2(n;3 - 1/2). 

Finally, the equations for the probability distribution of the number of photons in the 
modes a and /3, ^(72^) = = Y.np Pn^,n,3 and p{np) = p„^ = Pn^,ng, follow: 



-p(n„ - 1) - CiUaPM 



l + 52 + (5/A)n„"^ " ^ ^ ""^ l + 52 + (5/^)(^^ + l)^ 

+ Ci{na + l)p(na + 1) - 2C^raQ^ra^/r~^^_„^ [2p(ra„,ra;3) -^{71^,71^ - 1) - J9(ra„ - I,?!;?)] 



+ 2C3 (nc, + 1) X] ^P^nl+i,ng [2p(«a + 1, «/?) - p(n„ + 1, U/j - 1) - np)] 



(20) 



p(ra^) = -C2npp{np) + C2(r2/3 + l)p(n^ + 1) 
- 2Clnp ^ n„/s:;'^^^„^ [2p(n«, n^?) - p(nc, - 1, rip) - p{na, rip - 1)] 

+ 2C|(n^ + l)^n«K-;, 

n,3+l [2P('^Q) + 1) ^('^Q l7''^/3 + l) P('T'a, ■'T'/?)] • (21) 

In the steady-state regime, piria) = Pinp) = and Eqs. (!20|) and (I2T1) can be reduced 
to two-term recurrence relations by using the detailed balance condition and assuming that 
F{ni,nj)p{ni,nj) ~ F(nj, nj)p{ni). Here and denote the average photon numbers 
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in the modes a and /3, respectively. For n^, n/? > 1, the resulting equations are 



A 



- + 2ClnpK~\-\ = 0, 

J 



Vl + 52 + (i?/A)n, 
- 1) X 2Cln^Kzl^^ = 0. 



(22) 



(23) 



The equations (fT8|) - fl23|) are the main result of this work. Supplemented by the nor- 
malization condition pin^) = Xln^ ~ they will allow us to analyze the photon 
statistics, the threshold, the photon number fluctuations and the linewidth of the open-cavity 
laser in the steady-state regime. 



A. Photon statistics 



Photon number distributions p{na) and p{ni3) can be readily obtained by solving Eqs. 
(l22l) and (1231) numerically. But even without any numerical solution, it is easy to convince 
oneself that the only solution of Eq. (1231) is p(n^) = for n/j > 1 [and hence p(n/3 = 0) = 1 by 
normalization]. Equation fl23|) does not contain any gain, only damping terms are present. 
Besides, the terms proportional to C3 cancel each other well above threshold, i.e. for Ua ^ 1, 
and therefore the steady-state solution vanishes, i.e. pi'^n = 0, with the exception that 
p'^Q = 1 involving the normalization condition. In contrast, Eq. ( 1221) does have a non-trivial 
solution p{na) > for tIq, > 1 and this solution depends on the pump rate r. We therefore 
expect that if the laser effect occurs in our system, we should look for it signatures in the 
behavior of the composite mode a. 

In the limit of weak pump r — 0, we may consider the linear approximation for the 
laser equations, i.e. -8 = 0, and the photon number distribution resulting from Eq. ( I22l) 
approaches the thermal distribution as we see in Fig. 1. The analytical solution of Eq. ( 122|) 
(with n/3 ~ 0) can be approximated by 

with Ci = Ci(l + 5^), which is similar to the standard result for the single-mode laser jsj. 
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FIG. 1: Steady-state photon statistics for the composite mode a (dotted curve) below (curve 1), 
at (curve 2), and above (curve 3) threshold. The dashed and solid lines show the thermal and the 
Poisson distributions corresponding to the same average photon numbers fia as curves 1 and 3, 
respectively. For this figure, we fixed gi/T = 0.05, 5'2/r = 0.07, 5/T = 3, 711 /F = 6, 722/F = 5, 
and 712/F = 5.5. 

In the limit of strong pump, r — )■ 00, saturation effects become important and Bfia/A ^ 
1 + 5^. The analytical solution of Eq. ( |22l) tends to 



p{na) ^p(O) 



{A/B)\{A^/BCi 
(n„ + I/S)! 



(25) 



where A = A{1 +5^), and p{0) can be determined from the normalization of p{na). We thus 
observe that the distribution of Ua changes qualitatively when the pump is increased and 
that its limiting forms ( l24l) and ( l25l) coincide with those for the single- mode laser 8[. This 
suggests that the laser transition occurs for the composite mode a in our two-mode model. 
This is illustrated by the distributions of the photon number Ua found by solving Eqs. (!22|l 
and fl23l) numerically that we show in Fig. [H We see that at low pump (below threshold), 
the distribution is close to the thermal one: p{na) = (1 — exp[— ^w/Zc^T]) expl—nahw / ksT], 
where the effective temperature T is determined by exp[—hw/kBT] = A/Ci. At strong pump 
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FIG. 2: The average number of photons in the composite mode a as a function of the pumping 
rate r. The off-diagonal elements of the symmetric matrix 7 for the four curves are 712 = 0, 4r, 
Sr and 16r, respectively. Other parameters are as in Fig. [TJ The inset shows the dependence of 
the threshold on 712/r. 

(above threshold), pijia) approaches the Poisson distribution: p{na) = exp[—na]na' / i^a-) ^ 
where is given by Eq. ( I26l) below. 

B. Average photon number 

The average photon number fia can be found from Eq. ( l22l) as = '^aP('T-Q)- Far 
above threshold, the distribution of is strongly peaked around and p{na + 1) — p{na), 
as well as Kn^+i^np - Kn^^fip+i - Kn^^np- Together with = this yields 



Hence, the threshold for the composite mode a is given by the condition A/Ci = 1 + 6"^. 

The dependence of the threshold pumping rate r on the off-diagonal element 712 of the 
damping matrix 7 is shown in the inset of Fig. [2l In contrast, the mode (3 does not have 
a threshold and the number of photons in it is always equal to zero. The full dependence 




(26) 
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FIG. 3: The second-order correlation function (7^^^(0) as a function of the pumping rate, for the 
same parameters as in Fig. [2l 

of Ua on the parameters of the problem can be obtained by solving Eq. (|22l) numerically. 
In Fig. [2] we show the dependence of fia on the pumping rate r for different values of the 
off-diagonal element 712 of the damping matrix 7. Figures [1] and |2] show that the composite 
mode a behaves as a lasing mode with a well-defined threshold. This is also highlighted 
by the formal equivalence of Eq^ fl26l) and the standard expression for the average photon 
number in a single-mode laser 



C. Photon number fluctuations 

A common way to characterize fluctuations of the photon number in a mode of the 
electromagnetic field is to compute the so-called Mandel parameter, 

^ — , , 

Q = 1. 27 

n 

Equations and (12^ readily allow us to compute this quantity for the composite mode 
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a analytically well below and far above threshold, respectively, as 

~ — , below threshold, (28) 

Ci-A 

Ci 

Qa — above threshold. (29) 

A-Ci 

Another quantity which is often used to characterize fluctuations of the photon num- 
ber in experiments (see, e.g., Ref. 2J]) is the second-order correlation function ^''•^^(O) = 
= Qa/n + 1. Its value ranges from ^'^^•'(O) = 2 for thermal light to 
(^'•^^(O) = 1 for coherent laser light. Hence, the dependence of g^'^\Q) on the pumping 
rate shown in Fig. [3] allows one to identify the laser transition quite clearly. 



D. Laser frequency and linewidth 

Information about the frequency and the linewidth of light emitted by the two-mode 
open-cavity laser can be extracted from the off-diagonal elements of the density matrix, Eq. 



mi) (8|. Defining Pn^,np{ki,h) = Pnc,n^;nc+fci,n^+fc2 5 ^e follow Ref. [23[ and use an ansatz 
PnQ,n^(^i, ^2) = —fJ'{ki,k2)pna,nis{ki,k2). When we insert this into Eq. (fTSl) . it follows that, 
up to the lowest non- vanishing order in kik2 (valid for k <^ n), 

kf f A/{n^ + l) + 2B CA klC2 



8 \l + S^ + {B/A){n^ + l + ki/2) + {B/AAykl n„y Sn^ 
iA5ki/2 



(30) 



1 + P + {B/A){n^ + 1 + k^/2) + {B/AAfkl' 
The linewidth of the laser emission corresponding to the composite mode a is given by 
the real part of /i(l, 0), 



2Do, = - 
4 



A/{n^ + 1) + 2B Ci 



(31) 



[1 + 6^ + {B/A){n^ + 3/2) + {B/AAy 

For B/A <^ 1 and 6 = the linewidth reduces to 2LL = {A + Cij/Afia. This formally 
coincides with the result of the standard laser theory 8|, except for the definition of Ci, 
which includes an additional term, oc 712, in our case of the open-cavity laser. In Fig. HI 
we plot the dependence of the linewidth of mode a on the pumping rate r for r at least 
25% above threshold. The dependence of the linewidth on the off-diagonal element of the 
damping matrix 7 is shown in the inset for r = 100, which is far above threshold for the 
range of 712 shown in the figure. 
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FIG. 4: Linewidth of the composite mode a as a function of the pumping rate r for the same 
parameters as in Fig. [2j The inset shows the dependence of the hnewidth on 712/r for r = 100. 

As in the usual single-mode laser with a cavity of high quality factor j3, Isj], the finite 
linewidth of our two-mode open-cavity laser is due to the spontaneous emission of the active 
atoms. Previous studies have shown that in an open cavity, the emission line is further 



broadened by a factor called the Petermann factor 



of Refs. 



25|, due to the nonorthogonality of 



the cavity modes 26l |. However, our analysis here is based on the quantization procedure 



ll| that relies on the expansion of the electromagnetic field inside and outside of 



the cavity in terms of orthogonal modes. It is therefore interesting to check if our model 
reproduces the large Petermann factor expected from previous studies |26|. We define 



K 



DM 



(32) 



where the linewidth Da is considered as a function of 712 that quantifies the openness of 
the system in our model. Already from Fig. |4]we see that increases with 712. This is 
confirmed by Fig. [5] where K is shown as a function of 712 far above threshold. A simple 
analytical expression for K can be obtained aX B/A and for 5 = 0: 

C^i(7i2)[A + Ci(7i2)] 



K 



Ci(0) [A-Ci(7i2)] [A + Ci(0)] 
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(33) 



1 
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7i2/r 

FIG. 5: The Petermann factor i^T as a function of 712/r for the same values of parameters as in 
Fig. [2] and large pumping rate r = 100. 

Because 6*1(712) oc 712, we find that the growth of K with 712 is roughly hnear, which 
is in agreement with Fig. [51 Thus, our model reproduces the increased Petermann factor 
in open laser systems, despite the orthogonality of the basis in which we quantized the 
electromagnetic field. 

Finally, the imaginary part of yu(l,0) yields the shift of the laser frequency with respect 
to u: 

A = ^ (34) 

This equation shows that the frequency shift depends on the average photon number 
which, in turn, is a function of the off-diagonal element of the damping matrix 7. 



IV. CONCLUSION 



We developed the quantum theory of a laser with an open cavity. The openness of the 
cavity is mathematically described by a nondiagonal damping matrix 7. Assuming that only 
two modes of the "cold" cavity are allowed to participate in the competition for gain, we 
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have shown that the modes strongly interact with each other and that a composite mode 
(denoted by a here) is built up. This composite mode shows all the properties of a typical 
laser mode: threshold behavior, photon statistics evolving from thermal to Poissonian as 
the pumping rate increases, augmented photon number fluctuations in the vicinity of the 
threshold, and linewidth narrowing. At the same time, the precise behavior of the composite 
mode a at given values of parameters depends explicitly on the off-diagonal element 712 of 
the damping matrix 7. More precisely, an increase of 712 rises the lasing threshold and 
broadens the laser emission line. 

An important result that cannot be obtained from a semi-classical theory and needs the 
quantum theory developed in this paper to be understood is the broadening of the emission 
linewidth due to the coupling of the cavity modes through a common bath (Fig. [S]). This 
broadening should be accessible experimentally: indeed, measurements of the laser linewidth 
in a laser with two coupled modes were already performed in, e.g., Ref. It is interesting 
to note that in that work the mechanism of coupling was different and led to narrowing of the 
emission linewidth due to the so-called correlated spontaneous emission in a system of two 
excited levels that provide inversion for lasing on transitions sharing a common lower level 

n 

(see Ref. |21| for the theoretical model). In contrast to this, the two- mode open-cavity laser 
considered in this paper exhibits widening of the emission linewidth. This would certainly 
limit applications of such a laser for high-precision measurements but might be beneficial 
for other ap plic ations where low-coherence light is required (like, e.g., optical coherence 
tomography [28|). 

One of the possible applications of our theory may lie in the field of random lasers. In this 
case, the damping matrix 7 should be treated as a random matrix and our results should 
be averaged over the statistical distribution of its elements. However, to obtain results that 
can be directly applied to random laser systems, one has to generalize our analysis to the 
multi-mode case because the number of active modes in random lasers is expected to be 
large [l3|. 
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